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Abstract. - In usual superconductivity (SC), the pairs have zero total momentum irrespective 
of their symmetry. Staggered SC would involve, instead, pairs with a finite commensurate total 
momentum, but such exotic states have never been proven to be realized in nature. Here we study 
for the first time the influence of particle-hole asymmetry on the competition of staggered SC with 
Charge Density Waves (CDW) in a ferromagnetic medium. We obtain unprecedented situations 
in which CDW and staggered SC coexist. We also obtain cases of a SC dome near the collapse 
of a CDW state as well as cascades of transitions that exhibit remarkable similarities with the 
pressure phase diagram in UGe2 suggesting that SC in this material may be staggered coexisting 
and competing with a CDW state. 



The field of unconventional superconductivity (SC) is 
an extraordinarily rich source of challenging problems for 
fundamental and applied physics. High-T c cuprates, heavy 
fermion materials, borocarbides and organic SC are exam- 
ples of unconventional SC. A multitude of unconventional 

\ SC states have been proposed but only few have been 
proven to be realized in real material systems. Singlet 

\ or triplet unconventional SC states considered as realized 
so far are characterized by a zero total momentum of the 
pairs indicating in fact that the superfluid density is ho- 
mogeneous in momentum space. In the present Letter we 
suggest that the surprising ferromagnetic (FM) SC state 
of UGe2 is instead staggered, in which case the pairs have a 
finite total momentum. We also prove that staggered SC 
may coexist with charge density wave (CDW) states ex- 
plaining features of the pressure phase diagram in UGe2 
and opening a new perspective for the discussion of other 
unconventional SC. 

The unexpected discovery of SC well inside the itiner- 
ant FM state of UGe2 under pressure [1] and subsequently 
in ZrZn2 [2] and URhGe [3] represent a fascinating chal- 
lenge for our understanding of SC [4-12]. The kind of 
unconventional SC involved as well as the complexity of 
the pressure phase diagrams in UGe2 remain a puzzle. In 
fact, below Tc(0) = 52 K an almost fully polarized FM 
state is observed. Applying hydrostatic pressure, FM is 
suddenly eliminated at p C 2 = 1.5 GPa. Around p c i = 1.1 



GPa, SC appears and at the optimum p* = 1.2 GPa one 
finds T c (p*) < 1 K. This is well in the FM state where 
the Curie temperature and FM moment are still signifi- 
cant with T c (p*) 30 K, i.e. 60% of the original T c (0). 
In addition to SC, presumably another phase is present 
inside the FM state below T*(0) = 30 K, and T*(p) also 
decreases with pressure until at p* it hits the optimum 
T c (p) of the SC dome [8, 13]. The nature of the T* phase 
is not clear, one possibility that we adopt here is to as- 
sociate it with a charge density wave (CDW) [1,8, 14, 15], 
a scenario supported also by LDA+U calculations [16]. 
This would explain the associated experimentally observed 
heat-capacity anomalies [17] and jumps in the magneti- 
zation on crossing the T* phase boundary [13]. Finally, 
recent NQR experiments suggest that the T* phase sur- 
vives even below the SC transition dividing the SC dome 
into two parts: a low pressure part where SC and the T* 
order coexist and a high pressure part where there is no 
signature of the T* order [18, 19]. 

We explore in this Letter the competition of staggered 
SC (i.e. zone boundary SC) with CDW in a strong fer- 
romagnetic background, a situation that has never been 
considered before. In fact, all previous theoretical inves- 
tigations of SC in UGe2 considered zero momentum (or 
zone center) SC. Note that the possibility of the relevance 
of staggered SC states for some heavy fermion compounds 
has first been considered in the past by D.L. Cox and 
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coworkers [20] in the context of a Ginzburg-Landau ap- 
proach. Moreover, the SC pairing of spinlcss (or single- 
spin) fermions has been introduced in [21] and discussed 
within a Ginzburg-Landau theory. Here we start from a 
mean field BCS-type Hamiltonian 

H = ^Ck4 c k - ^(WkC^Ck+Q + h.c.) 
k k 

-£(A£ctcL k + *.c.) -E( A ?4ci k _ Q + M (i) 

k k 

The first term describes a 2D tight binding FS whose nest- 
ing properties with a wave vector Q=(-7T, 7r) are controlled 
by the ratio of n.n. (ti) and n.n.n. (t2) hopping matrix el- 
ements. For t 2 /ti < 1 the nesting with Q is removed com- 
pletely and the FS changes its shape. This is a schematic 
model for the destruction of nesting of the LDA+U Fermi 
surface [16] under pressure. There are several possibili- 
ties of competing SC and CDW orders described by the 
gap functions W and A in Eq.(l). They are related to 
the effective interactions via gap equations like Eqs. (2,3). 
The effective interactions Vf c k , and V^jB^ of the itinerant 
5f-quasiparticles have a purely electronic origin. We may 
have both unconventional SC with zero total pair momen- 
tum A k and at finite pair momentum A k . The CDW gap 
function is denoted by Wk and like A° or A k belongs to 
an irreducible representation of the tetragonal D^ group 
(this is also the approximate symmetry of UGe2). 

To treat both SC and CDW order parameters in a com- 
pact manner we introduce a Nambu-type representation 
using the spinors ^ = (c k , c_k, C k+Q) c_k-cj) ■ Accord- 
ingly we use the tensor products p — (a ® I) and a = 
(7 (g) a) for the Nambu representation of the Hamilto- 
nian in Eq. (1). We assume that nesting in the fully 
FM polarized band is responsible for the CDW transition 
associated with the T* line in UGe2. Pressure reduces 
T* because it relaxes the nesting conditions. To model 
this effect we write the electron dispersion as a sum of 
particle-hole symmetric terms responsible for nesting and 
particle-hole asymmetric terms that represent the devia- 
tions from nesting: Ck = 7k + ^k where 27 k = £ k - £k+q 
and 25k = £k + Ck+Q • When Sk — there is particle-hole 
symmetry or perfect nesting with wavevector Q. Applica- 
tion of pressure adds a 5^ term in the dispersion in addi- 
tion to the 7k term already present at zero pressure. We 
classify the SC and CDW order parameters with respect 
to their behavior under inversion (I) k — > — k, translation 
(tq) k — > k + Q and time reversal (T) in the charge sector. 
Instead of the latter we may also use complex conjugation 
(C) which satisfies the equivalence relations C = -T (A k ); 
C ee IT (A k ) or C ee t Q (W k ). These discrete transfor- 
mations may then be used to classify the possible groups 
of competing SC/CDW order parameters. Obviously C is 
redundant for the three order parameters considered, but 
we include it in the notation for clarity. 

Because the spins are frozen, the q = SC pair states 



translation we have both signs A° + q = ±A° and un- 
der C we get (A°)* = - (A°) T = - A°. SC pair states 
with finite momentum may in principle have both par- 
ities: A Q „ = ±A k because the required antisymmetry 
may also come from the shift by a lattice vector R with 
cxp(iQR) = — 1. On the other hand the tq transla- 
tion requires that always A k+ q = — A k and under C we 

have (A^)* = (A Q k ) T = - A Q k . These transformation 
properties allow four possible SC order parameters, two 
at zone center and two at zone boundary or staggered SC: 



,01- 
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^k ! 
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A k Q7+ - 



Here the first in- 



dex or Q indicates the total momentum of the pair, the 
second index R or I indicates whether the order parame- 
ter is real or imaginary, the third index ± indicates parity 
under inversion I and the last index denotes gap symmetry 
under tq. As mentioned the index I(R) is redundant. 

For the CDW order parameter both odd and even states 
under I and tq arc allowed so that W_k = ±Wk and 
Wk+q = ±Wk may hold. Since C ee tq for this or- 
der parameter the redundant index R or I is associated 
with the tq-indcx ± respectively As a result we have 
here again four different possible (CDW) order param- 
eters: W R++ , W£ + ~, W£- + , W k ~~, where the in- 
dices have the same meaning as the last three indices 
in the SC order parameters. According to the above 
symmetry classification there are sixteen possible pairs 
of such competing SC/CDW states and only eight of 
them concern staggered SC on which we are interested 
here. Within our formalism we can calculate Green's func- 
tions and self-consistent gap equations for each of these 
eight cases. As an example relevant for UGc 2 we re- 
port here for the case of the competition of W k + ~ with 
A k R where the Hamiltonian in spinor representation 
with Pauli matrices di and % is H = J2k ^'Ipk^'L where 



Sk = -fkP3<5 3 + 6 k <7 3 + A^ R p 2 o- 2 
that if instead of having the competition of W^ + ~p2 with 
A k R P2&2 as above, we had any of the other pairs of 
competing order parameters, we would just have to replace 
the corresponding SC and CDW terms in the above Hamil- 
tonian. The Green's functions that result would be modi- 
fied accordingly. For the W£ + ~{>2 and A k R P2&2 order 
parameters, the most obvious realization in D 4 /j symme- 
try is a d-wave CDW and p-wave (finite momentum) SC 
order parameter given by A k = A^(sinfc x + smk y ) (i.e. 
E u (l, 1)) and Wk = W (cos k x cos ky) (i.e. Big) ■ From 
the Hamiltonians we obtain Green's functions and then 
self-consistent gap equations for both order parameters 
which after analytic summation over the Matsubara fre- 
quencies take the form of the following system of coupled 
equations that are reported for the first time here: 



W k + ~p 2 - We note 
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may only have odd parity with A°, 
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v k,k' 



hA k , 



tanh „ ' - tanh 



2T 



, £+( k ') 
tanh — ±^ — - + tanh 

2T 



2T 
£_(k') 



2T 



(2) 
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25±(k) 



(4) 



Here the effective potentials Vj? ^, , Vj? jB w are separable for 
the asummed E u and Bi g channels. If solutions of these 
coupled SC/CDW gap equations exist they are unique. 
For uniform order parameters assumed here they also have 
lower free energy as compared to the normal state [22,23]. 
Eqs. (2,3) account for the following four pairs of compet- 
ing CDW and zone boundary states: with W R++ , 
A R Q- with W I+ ~, A IQ - with W R++ and A 1 ®- with 
W I+ " . There is a second system of coupled gap equa- 
tions that describes the competition of the remaining four 
pairs of SC and CDW gaps: A R Q~ with W R ~+, A R ®- 
with W 1 — , A /( 2- with W R -+ and A 1 ®- with W 1 " : 



+ - 



AE + (k')B(k) 
S(k)- 7 ^ 



tanh 



2T 
E_(k>) 



4E_(k')B(k) 
^ = E^ k '{J^|ytanh 

fl(k)-5£, 



+ 



4E_(k')B(k) 



tanh 



2T 

2T 
E_(k>) 



2T 



(5) 



(6) 



where 



E ± (k) = 



AlWl_ 

wl 



+ 




(7) 



S(k) = ^(^+W2,)+7M, (8) 



We have solved selfconsistently the systems of Eqs. (2,3) 
and (5, 6) for a 2D tight-binding model on a square lattice. 
In that case, the particle-hole symmetric term corresponds 
to nearest neighbor hoping 7k = ti(cosfc x + cosk y ) while 
particle-hole asymmetry is introduced by the next-nearest 
neighbor hopping terms 5^ — £2 cos k x cos k y . We have 
performed a large number of self consistent calculations 
varying the pairing potentials in the two channels pro- 
ducing eight maps (two of them reported in figure 1) of 
all possible transitions induced by particle-hole asymmetry 
(i.e. by pressure) in the low-T region for all the pairs of 
competing CDW and staggered SC order parameters that 
are possible. To take into consideration the fact that the 
different CDW and SC gap symmetries involved may cor- 
respond to different momentum structures for the order 
parameters, we have considered the separable potentials 
approximation that allows to search for solutions of a spe- 
cific momentum structure. Therefore, the axes in Fig.l are 
the amplitudes y CDW and V sc of the pairing interactions 
and it is understood that a corresponding form factor has 
been considered. We have investigated the coexistence of 
order parameters in the whole (y CDW , V 5C ) plane from 
the moderate coupling to the strong coupling regime since 
we have no microscopic derivation for the effective pairing 
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Fig. 1: Maps of the dependence of phase sequences on the ef- 
fective interactions y CDW and V sc for low temperature. Ar- 
rows indicate the cascade of phases obtained when £2 /ti grows 
starting from zero. The black dots separate regions of differ- 
ent phase sequences under growing tijt\. All phases coexist 
with ferromagnetism (FM). The phases indicated as FM, are 
phases in which there is not any finite SC or CDW order pa- 
rameter and so only FM is present. Figure a) corresponds to 
the competition of A R ®~~ with W 1 . Figure b) corresponds 
to the competition of A R ® with W R ~ + . The potentials are 
in units of ti. 



strengths. Arrows in Fig. 1 indicate the cascade of phases 
observed when the ratio £2/^1 grows starting from zero. 
Since we consider a spin polarized background, all states 
reported also coexist with FM, and the transitions to the 
FM state reported at high values of r^Ai has the meaning 
of a transition to a state that is only ferromagnetic with 
no CDW or SC order parameter present. We note in figure 



la that in a large portion of the V , V 



SC \rCDW 



parameter 



space there is at low temperature a transition from a CDW 
state to a state in which CDW and SC coexist. The same 
transition is also present over a portion of the parameter 
space in the case of figure lb. We note that the coexis- 
tence of zone-center SC with CDW has been reported in 
previous theoretical studies [22,27,28]. 

We now look more closely to the situations in which 
CDW and staggered SC may coexist. We report in figure 
2a the behavior of the CDW+SC state with tijt\ in a 
characteristic example that corresponds to V CDW = 4 
and V sc = 2.5 in the competition of A RQ with W R ~ + 
at low-T the mapping of which is reported in figure lb. We 
observe the cascade of transitions from CDW to CDW+SC 
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Fig. 2: (Color online): A characteristic example in which devi- 
ations from nesting induced by fe/ii lead to coexistence of the 
W R ~ + CDW (open circles in red) with the A RQ ~~ staggered 
SC order parameter (filed circles in black) corresponding to 
yCDw _ 4 anc j y-sc _ 2 5 ODServe the cascade of transi- 
tions from CDW to CDW+SC and then to FM as was already 
reported in figure lb in this region of of y CDW and V sc val- 
ues. In the inset is shown the temperature dependence of both 
order parameters when tz/ti = 0.9 in which case the CDW and 
SC orders coexist at low temperatures. All quantities are in 
units of ti. 



at t 2 /h w 0.72 and finally to the FM state for t 2 /h > 0.93 
in agreement with the cascade reported in figure lb for 
these couplings. It is remarkable in figure 2a that the 
transition from CDW to CDW+SC is smooth as a function 
of ^Ai in the low-T regime whereas the transition from 
the CDW+SC state to the FM state (i.e. the state with no 
SC or CDW) is first order in t-xft\. In figure 2b we show 
the corresponding transitions with temperature when we 
take tijt\ — 0.9. In this case, according to figure 2a, 
we have indeed in low-T coexistence of CDW and SC. 
We observe the counter-intuitive behavior that when SC 
appears as we lower the temperature, the CDW gap grows 
instead of being reduced. 

A particularly interesting cascade of transitions in rela- 
tion to the observations in UGe2, is the one from CDW 
to CDW+SC to SC and finally to FM. This cascade is 
observed over a large portion of the parameter space of 
pairing potentials when the A R ® and W compete 
(cf. fig. la) and over a smaller portion in the case of 
competition of A R ® with W R h (fig. lb). A cascade 
of transitions in the low temperature regime that exhibits 
amazing similarities with that observed in UGe2 as a func- 
tion of pressure is shown in Fig. 3. It corresponds to the 
competition of W R h CDW with the staggered SC order 
of the form A RQ — wne n V sc = 3.5 and V CDW = 12.5. 
The maximal critical temperature of SC coincides with 
the crossing of the CDW critical line. Moreover, the SC 
dome is divided into a part in which SC and CDW coex- 
ist and a part in which only SC is present (with FM of 
course). Quite remarkably the SC critical temperature is 
reduced almost linearly at the higher values of t 2 /t\. The 
above results are in surprising qualitative agreement with 



Fig. 3: Evolution with temperature of a cascade of t2/t\ in- 
duced transitions from CDW to CDW+SC to SC and then to 
only FM obtained when V CDW = 12.5 and V sc = 3.5. The 
SC and CDW order parameters considered are A RQ ~- and 
W R This phase diagram shows striking similarities with 
features of the pressure phase diagram of UGe2 if one suppose 
that the T x phase corresponds to a W R CDW ordering and 
of course SC to A RQ — . 

findings in UGe2 [18]. In particular, recent NQR results 
indicate that the T* phase coexists indeed with SC over a 
portion of the SC dome [19] and our results are the first 
to provide a theoretical picture for it. 

Staggered SC states are relevant for magnetic SC [20,24, 
25] because they are similar to the Fulde-Ferrel states [26] 
except that the modulation of the superfluid density co- 
incides with the characteristic wavevector of the CDW. It 
appears, therefore, plausible to consider these states in the 
analysis of SC in UGe2 which is observed only in the FM 
regime because for a staggered SC the FM background is 
necessary in the same way as the magnetic field is neces- 
sary in order to obtain the usual Fulde-Ferrel phase. Since 
pressure eliminates the FM state it naturally eliminates si- 
multaneously the staggered SC state as well. Note finally 
that we obtain staggered SC states only over a limited 
dome near the collapse of the CDW phase as in figures 2 
and 3 within a mean field approach without any fluctua- 
tions involved. 

In conclusion, we have demonstrated the possibility to 
have coexistence of staggered SC with CDW and cascades 
of transitions induced by particle-hole asymmetry that re- 
produce the pressure phase diagram observed in UGe2 
identifying the T* phase as a CDW phase. Such exotic 
states may be relevant for other magnetic SC as well. 

* * * 
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